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Ch . Abstract 

Low-energy El strength in neutron-rich "'^'^^^^^'^Sn isotopes is analyzed in the framework of the Skyrme 
^ ' random phase approximation (RPA) with different Skyrme forces. A double folding procedure is applied 

i> : 

Q>^ , to take into account the energy-dependent width effects beyond RPA. All the considered Skyrme forces 



indicate a soft prolate deformation in the open shell isotopes ^'^^^^^"'Sn. The integrated El strength in the 
energy region of the pygmy resonance grows with the neutron number. The influence of deformation on 
the integrated strength near the particle emission thresholds (which is of a keen astrophysical interest) is 
strictly suppressed by the mutual compensation effect for the branches of the giant dipole resonance. The 
results obtained are in a good agreement with the previous findings of the relativistic mean field model. 
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I. INTRODUCTION 



Exotic nuclear collective excitations, like the dipole pygmy mode (PM), represent a subject of 
intense investigations during last decades, see the recent review [l^. The analysis of these excitations 
provides a deeper insight into nuclear dynamics. Besides, the PM affects the El strength near the 
particle emission thresholds, which is of a key interest for astrophysical applications The 
.-•M is viewed as a vibration of the neutron skin against the isospin-saturated proton-neutron core 
1|. In this connection, the evolution of PM with rising the neutron number and approaching the 
neutron drip line attracts a large attention. 

The available experimental data show the growth of PM strength with neutron excess in 



medium- heavy Q-5] and light [s 



, Ml stable nuclei. This trend was confirmed by relativistic and 

n 

non-relativistic (Skyrme) mean-field calculations Meanwhile, the experimental effort focuses on 
drip-line regions and thus a theoretical analysis of the PM behaviour in unstable nuclei is required. 
Such study was recently performed for neutron-rich ^^^~^^^Sn isotopes within the relativistic mean- 
field model (RMF) P|. A growth of PM strength with the neutron number was confirmed and a 
region of deformed isotopes ^'^^~^^^Sn with more spread of the strength was predicted. 

The aim of present paper is to investigate the deformation and neutron-excess effects in 
the ^^^~^^^Sn chain within the Skyrme mean-field approach. The separable random-phase- 
approximation (SRPA) model 0, [lO] is used. SRPA is fully self-consistent as both the static 
mean field and residual interaction are derived from the same Skyrme functional. The model cov- 
ers both spherical and axially-deformed nuclei. It was already applied to describe electric jgl 1^ 
and magnetic [isl. \l\ giant resonances as well as the El strength near the particle thresholds 15]. 



In the present paper, the SRPA formalism is supplemented by the double folding procedure which 
allows to compute strength functions with an energy-dependent width. The SRPA analysis below 
is done with the representative set of Skyrme forces, SkT6 Q, SkM* 0, SLy6 Q, and SkI3 
[3], covering a broad interval of Skyrme parameters and nuclear matter characteristics. 



II. SRPA APPROACH 



SRPA is derived from the Skyrme energy functional 

E{p,T,s,3,T,x) = I n{r)dr (1) 

involving time-even (nucleon pq, kinetic-energy Tg, spin-orbit Jq) and time-odd (current jq, spin 
Sq, vector kinetic-energy T^) densities, as well as the pairing density Xg- The index q stands for 



neutrons {q = n) and protons {q = p). 

The energy density 'H(r) consists of the kinetic-energy ^kin(r), Skyrme 'Hsk(r), pairing 'Hpair(r), 
and Coulomb ?^coui(r) contributions 



n{v) = -HMnir) + nsk{r) + ?^pair(r) + ^Coul(r 
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where the densities without index, hke p = pp + pn, are total. The explicit expressions for the 



densities can be found elsewhere 



10( 1 . In the part ([5]) , is the pairing strength and po is the equi- 



librium nuclear matter density. The pairing is treated at the BCS level. The Coulomb contribution 
includes direct and exchange terms. In the part (jH), the Skyrme terms with the tilded parameters 
hi and 6- are relevant only for magnetic excitations and so are omitted in the present study. 
The self-consistent derivation P, |l^ gives from (HJ-O the SRPA Hamiltonian 



H = /iHFB + Vra 



(7) 



where /ihfb is the HFB mean field 



^HFB = j dry^ 



dE 



dJa+ (r) 



(8) 



and V^cs is the separable residual interaction 

1 ^ 

2 ^2 I'^kk'^kXk' + rjkkiYkYki 



k,k'=l 



(9) 



with one-body operators 



Xk = i j dr j dr' 
Yk = i j dr j dr' 



dJa+ir')dJa+ir) 
9J„_(r')5J„_(r)J 



PkJaM |)J«'(r'), 



Qfc,J,_(r) |)J«'_(r') 



(10) 



and inverse strength matrices 



^kk' 



(11) 



Here q+ = p, r, J,x and a_ = j,s, T enumerate time-even Ja+ and time-odd Ja_ densities; Ja± are 
the density operators; and P ^; = i\H , Qk] are time-even and time-odd hermitian input operators 



chosen fohowing the procedure 



Id, 



12]. The operators of the residual interaction and Yk are 
time-even and time-odd, respectively. The single-particle Hamiltonian /ihfb is determined by the 
first functional derivatives of the initial functional ([1]) while operators X^ and are driven by the 
second functional derivatives of the same functional. Hence the model is fully self-consistent. The 
number K of separable terms is determined by the number of the input operators Qk- Usually we 
have K = 3 — 5. Hence a low rank of the RPA matrix and minor computational effort even for 
heavy deformed nuclei. 

When studying giant resonances, one is usually interested in the strength function. For electric 
dipole excitations, the strength function reads 



Sei{E)=J2Sei^.{E), Sei^{E) = Y,EumfEif.\)\^CA{E - E, 



(12) 



where /eiii = J2i ef'^erYi^ is operator of £'l/i-transition with the effective charges = N/A and 



eff 



„efT 



-Z/A and 



(AiE-E, 



1 



A 



2^ (^E-E^f + i^y 

is the Lorentz weight with the averaging parameter A to simulate the smoothing effects beyond 
SRPA (escape widths and coupling to complex configurations). The strength function is composed 
by computing the detailed RPA spectrum E^, and reduced transition probabilities B{El,v) = 
K^I/bi^DP- However it is more convenient to compute the strength function directly by using 
the linear response techniques [3] for which the factorization Q additionally simplifies the 
calculation. We thus use this direct technique. 

The actual smoothing effects depend on the excitation energy E. The escape widths appear 
above the particle emission thresholds and grow with energy due to widening of the emission phase 



(13) 



space. The coupling with complex configurations (collisional width) also increases with energy. 
To take into account these trends, one should use in P^ - (fT3|) an energy-dependent averaging 
parameter A{E). This can be done by implementation of the double folding scheme which allows 
to use the strength function formalism (with fixed width) for the RPA solutions. At the first step, 
we calculate the strength function (jl2|) with a small but fixed value of A. This gives the strength 
distribution Sei^i{E') very close to the actual RPA one but for the equidistant energy grid. At the 
next step, this strength is additionally folded by using the energy dependent averaging A(i?): 

S'e^{E) = j dE' Sei{E')Ca(E'){E - E') . (14) 

In the present study we use the simple linear energy dependence 

f Ao for E' < Eth 

A{E') = { - (15) 

[ Ao + a{E' - ^th) for E' > Ett 

where i?th is the energy of the first emission threshold and Ao is a minimal width describing the 
coupling with complex configurations below the threshold. The rate a is chosen to give a proper 
reproduction of the form of the giant dipole resonance (GDR). In the present study, we use Ao = 0.1 
MeV and a = 1/6. 

The dipole photoabsorption cross section is related to the strength function ()12p as cr{E) ~ 
0A02SEiiE) with SsiiE) being in e^fm^ and a{E) in fm^. 

III. RESULTS AND DISCUSSION 



The results are obtained for the representative set of Skyrme forces: SkT6 flffl, SkM* 171, 

n □ t J, ujj, 

SLy6 [18[, and SkI3 [19|]. These forces have different nuclear matter characteristics and are widely 
used for the description of ground state properties and dynamics of atomic nuclei [2^, including 
deformed ones. The calculations use a cylindrical coordinate-space grid with the mesh size 0.7 fm. 
The equilibrium quadrupole deformations are found by minimization of the total energy of the 
system. 

In Fig.l the calculated equilibrium quadrupole deformations /3 in Sn isotopes with 100 < A < 
166 are shown for the four different Skyrme forces. The broad mass region involves neutron- 
deficient, stable, and neutron-rich nuclei. The deformations are given for the total, neutron, and 
proton densities. The results are compared with the RMF findings jsl. It is seen that for A < 140 
the isotopes keep the spherical shape, for exception of a narrow region 112 < A < 116 where forces 
SkT6, SLy6, and SkI3 predict some oblate softness. It is remarkable that, despite the semi-magic 
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FIG. 1: Equilibrium total, neutron and proton quadrupole deformations in ^^'^ ^^^Sn, calculated with 
Skyrme forces SLy6, SkM*, SkT6, and SkI3 and compared to the RMF resuhs 




FIG. 2: Dependence of the HFB nuclear energy on the total quadrupole deformation /3 for SkT6, SLy6, 
SkM*, and SkI3 forces in selected neutron-rich Sn isotopes. 

character of Sn isotopes (Z=50), they exhibit a significant prolate deformation in the neutron- 
rich region 142 < ^1 < 162, which is in accordance to RMF predictions [8]. Note that neutron 
deformation in this region is systematically larger than the proton ones, which may signify the 
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FIG. 3: The doubly folded El strength function for SkM* and SLy6 in selected Sn isotopes. 

higher deformation of the excess neutrons. Perhaps, just the drive of the large neutron skin to the 
deformation causes the whole effect in these semi-magic nuclei which otherwise might be spherical 
as in the neutron-deficit region. 

The further analysis shows that ^^^~^^^Sn isotopes should be rather considered as prolate-soft 
instead of prolate-deformed. Indeed, the total HFB energies ii^HFB depicted in Fig. 2 demonstrate 
for 142 < A < 162 very shallow deformation minima. Though the deformations and -Ehfb values in 
the transition regions with A=142-146 and 158-162 noticeably vary with the force (in some cases, e. 
g. for SkI3 in ^^^Sn, both prolate and oblate minima are predicted), the results for different forces 
are qualitatively similar. Altogether they point out to possible strong softness and anharmonicity 
in 142-I62sn. 

In Fig. 3, the doubly folded strength function S'^i{E) is shown for the forces SkM* and SLy6 
in the representative isotopes (spherical ^^^Sn, transition ^'^^Sn, deformed ^^^Sn, and spherical 
^^^Sn) covering the deformation region and its vicinities. The double folding is computed with the 
threshold energies (associated to the neutron single-particle Fermi levels) marked in Fig. 4. It gives 
a reasonable width A ~ 1 MeV at the GDR peak. Due to larger averaging at higher energies, the 
double folding considerably suppress the artificial strong shoulder at the right GDR side, which 



often appears in Skyrme-RPA calculations both with [S 



11 



12j | and without {2I] the separable 
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FIG. 4: The low-energy El strength function (with constant A = 0.1 MeV) for SkM* and SLy6 in selected 
Sn isotopes. The estimated particle emission thresholds Eth are marked by arrows. 



prescription. Fig. 3 shows that SkM* and SLy6 give similar shapes for spherical and deformed 
isotopes but rather different deformations in transition ^^^Sn. Some fine structures are seen in the 
pygmy energy region 5-10 MeV. The El strength in this region grows with the neutron number. 

The latter effect is better visible in Fig. 4 where the El strength function (jl2p is exhibited 
with a small constant averaging A = 0.1 MeV in the low-energy pygmy region. It is seen that 
fine structure of the strength significantly depends on the force and it is more separated from 
the GDR in spherical isotopes than in deformed ones. Hence, in accordance to RMF results [s], 
spherical nuclei are indeed more suitable for observation of the pygmy mode. Note also that pygmy 
energy regions stay near the particle emission thresholds and so the pygmy resonance may affect 
the near-threshold El strength of the astrophysical interest 

In Fig. 5, the pygmy nature of the low-energy strength is illustrated for particular RPA states 
in spherical ^^■^Sn and ^^^Sn. It is seen that in the nuclear surface at 6-8 fm the neutron transition 
density dominates thus manifesting the oscillation of the neutron excess against the proton-neutron 
core. 
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FIG. 5: Neutron and proton transition densities Pv(r) in particular low-energy RPA j/-states with large 
B(E\^v) values in spherical -'^'^^Sn and ^^''Sn. The force SLy6 is used. 

Finally Fig. 6 shows the integrated El photoabsorption cross section 

r-E 



a{E')dE' 



4MeV 



(16) 



in units of the Thomas-Reiche-Kuhn (TRK) sum rule YItrk ~ 60A^Z/A mb MeV. The upper 
panel compares the ratios in spherical ^^^Sn and deformed ^^^Sn. At E > 7 MeV the isotope 
"^^^Sn delivers more low-energy El strength. This may be caused by two factors in ^^^Sn, the 
deformation and the larger mass number. The deformation down-shifts the left ^=0 GDR branch 
and the larger nuclear size down-shifts the GDR in general. While the latter factor is clear, the 
effect of deformation deserves more analysis {l5 . 



22 



23l ]. It is inspected for ^^^Sn in the lower 



panel of Fig. 6. The comparison of the photoabsorption cross sections for equilibrium deformation 
/3 = 0.316 and deliberately constrained spherical case /3 = shows that the deformation results in 
a larger cumulative strength only at > 10 MeV, i.e. in a close vicinity to the GDR. Instead, at 
the lower energies near thresholds (the region of the astrophysical interest) the deformation leads, 
in fact, to a weak depletion of the strength. 

The weak deformation effect in the pygmy region can be explained by the destructive competi- 
tion of /i = and ^ = 1 GDR contributions to the low-energy El strength [3]. As seen from Fig. 
6, both these GDR branches affect the pygmy region and the /i = branch evidently dominates. 
However, at low energies far enough from the GDR (E < 10 MeV in our case), the influences of 
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FIG. 6: The integrated photoabsorption cross section ([TE)) in the percentage of the TRK sum rule Strk- 
The upper panel compares the cross sections for spherical ^"^^Sn and deformed ^'^'*Sn. The bottom panel 
shows the cross sections in -'^^^Sn for the equilibrium deformation /3 — 0.316 and constrained spherical case 
/3 = 0. The contributions of /i = and /i = 1 GDR branches for /3 = 0.316 are exhibited by areas with linear 
and check filling, respectively. The force SLy6 is used. 



both branches become more comparable. Indeed, the fi = branch, being closer to the pygmy 
region, carries at the same time twice smaller strength than the = 1 one. Hence the comparable 
effects of the branches at the remote energies. The deformation shifts the ^ = and = 1 branches 
in the opposite directions, thus leading to the deformation contributions of the opposite sign and 
comparable magnitude. These contributions compensate each other and drastically suppress the 
total deformation effect at the remote low energies. In principle, the rest effect can be of any 
sign, e.g. negative, as in our case. A weak negative deformation effect agrees with findings [8] and 



contradicts the assertion 



23j] on its strong and positive magnitude. 



IV. CONCLUSIONS 



The long chain of neutron-rich Sn isotopes with 132 < A < 166 is investigated within the 
Skyrme RPA method by using SkT6, SkM*, SLy6, and SkI3 Skyrme forces. Both spherical and 
axially deformed shapes are considered. The smoothing effects beyond RPA are simulated by 
the double folding procedure which allows a direct computation of the strength function with the 
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energy dependent width. 

All the Skyrme forces predict a region of prolate-soft isotopes ^^^~^^^Sn for which a strong 
anharmonicity in low-energy collective modes is expected. Though prolate shapes look surprising 
for the semi-magic Sn isotopes, they can be caused by the large neutron excess driving the system 
to deformations. The calculations also show that the pygmy El strength increases with the neutron 
number. Following our analysis, the nuclear deformation does not noticeably affect this trend and 
leads mainly to more spread and redistribution of the El strength. The effect of deformation on 
the cumulative El strength near the particle emission thresholds is shown to be small and negative 
(slight depletion of the strength) . The result is explained by the compensation of the deformation 
contributions from = and = 1 GDR branches. This finding may be of a significant interest 
for astrophysical applications. 

The four Skyrme forces yield qualitatively similar results though their numerical predictions 
can vary for the total nuclear energy and fine structure of the pygmy mode. Our results are in a 
good agreement with the RMF study [8] 
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